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ABSTRACT 

A gamma ray burst with strong optical-UV emission occuring in a molecular cloud 
will photodissociate H2, photoionize H2, H, and He, and destroy dust grains. We model 
these processes, including time-dependent radiative transfer in both continuum radia- 
tion and the resonance lines of H2. The UV will pump H2 into vibrationally-excited 
levels. We calculate the absorption spectrum imprinted on radiation from the GRB at 
various times. In addition to the strong absorption lines of f = H2 at A < 1110 A due 
to cold ambient gas, we find that radiation reaching us from the GRB and its afterglow 
will show strong absorption lines due to vibrationally-excited H2 at 1110 < A < 1705 A. 
These absorption lines, if observed, would provide unequivocal evidence for associa- 
tion of the GRB with molecular gas. Low-resolution spectra will exhibit conspicuous 
features due to clustering of individual lines; a list of the strongest such absorption 
features is given for spectral resolution R ~ 350 characteristic of the grism on the Swift 
UV-Optical Telescope. 

Subject headings: galaxies: ISM ~ gamma rays: bursts - molecular processes - ISM: 
clouds - ISM: molecules 

1. Introduction 

Gamma-ray bursts (GRBs) rank among the most dramatic, most energetic, and least un- 
derstood celestial phenomena. Although at least some GRBs are now known to originate from 
cosmological distances, the progenitor objects, or their locations in galaxies, are not yet known. 
Only ~40% of GRBs have detectable optical afterglows, whereas the X-ray afterglow detection rate 
is ~100%; the undetected optical afterglows may have been extinguished by dust in the host galaxy 
(see Lazzati, Covino, & Ghisselini 2001; Ramirez-Ruiz, Trentham, & Blain 2001; and references 
therein). Since GRBs may be associated with star- forming regions (Paczyiiski 1998, 1999; Mac- 
Fadyen, Woosley, & Heger 2001), it is of interest to consider phenomena which would be indicative 
of molecular gas in the vicinity of the GRB. 
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At least some GRBs are accompanied by intense optical emission, as demonstrated by detection 
of a 9th magnitude optical transient coinciding with GRB 990123 (Akerlof et al. 1999) and optical 
afterglows associated with other GRBs (e.g., GRB 990510: Stanek et al. 1999; Israel et al. 1999). 
The energy radiated in the optical- UV flash and afterglow can be substantial, and can have dramatic 
effects on interstellar gas and dust in the vicinity of the GRB. The hv > 13.6 eV emission will 
photoionize the gas, and the pulse of optical radiation will vaporize dust grains out to substantial 
distances from the GRB (Waxman & Draine 2000, hereafter WDOO). The gamma-rays and hard 
X-rays emitted by the GRB will contribute to ionization of the nearby gas, but the UV and soft 
X-rays have the dominant effect because of the much greater number of photons, and much larger 
photoabsorption and photoionization cross sections. The time-dependence of atomic and ionic 
absorption lines and photoionization edges in the spectra of GRBs has been discussed by Perna & 
Loeb (1998), Bottcher et al (1999), and Lazzati, Perna, &; Loeb (2001). In the present work we 
focus on the rich absorption spectrum of H2 . 

Ultraviolet radiation will destroy H2 in the vicinity of the GRB, but before destruction some 
of the H2 molecules will be vibrationally excited by UV pumping. A preliminary study of the UV 
pumping (Draine 2000) found that high degrees of vibrational excitation would prevail in the II2 
undergoing photodissociation and photoionization, and an estimate of the column density of this 
vibrationally-excited gas indicated that it could imprint a conspicuous absorption spectrum on light 
reaching us from the GRB flash and afterglow. If observed, this absorption at A < 1650 A would 
be a clear sign that the GRB was in close proximity to a molecular cloud. 

The previous study did not explicitly include self-shielding in the H2 lines, or the details of 
the competition between photoionization and photodissociation of the H2, and therefore it was 
only possible to make approximate estimates of the column density of vibrationally-excited H2. 
In the present work we carry out a detailed calculation of the effects of a high-luminosity optical 
transient on surrounding dust and molecular gas. We calculate the evolution of the dust and gas 
in the neighborhood of the source of radiation, inchiding radiative transfer through the absorbing 
medium. We include the dominant sources of opacity: dust grains, H2, H+, H, He, and He+. 
Computational limitations dictate using radial zones which can be optically thick to 13.6 eV - 
100 eV ionizing photons, and in the Lyman and Werner band transitions of H2. We describe a 
numerical method which can be used to follow the motion of the ionization and dissociation fronts 
in this regime. 

The GRB light curve adopted for the calculation is given in §2, where it is compared with the 
observed light curve of GRB 990123. Because GRB 990123 appears to have been unusually bright, 
in this paper we consider the effects of an optical transient with a peak ultraviolet luminosity less 
than 1/40 that of GRB 990123, and an integrated optical luminosity per steradian perhaps 1% 
of GRB 990123. The computational method is described in §3, and the physical processes (dust 
destruction, photoexcitation, and photoionization) are summarized in §4. 

In §6 we present results for our adopted lightcurve Li,{t), for two cloud densities (nn = 10^ and 
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cm~^) and for two different power-law indices {P = —0.5 and -1) for the UV- X-ray spectrum 
Ljy oc u^. In §7 we discuss the prospects for observing these absorption features in the spectrum of 
a GRB firebah or afterglow, either with ground-based observations of GRBs at redshift z > 1.25, 
or using the grism on the Swift GRB Explorer for GRBs at redshift z > 0.06. We summarize our 
results in S8. 
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Fig. 1. — Filled symbols: observed light curve of GRB 990123 (Akerlof et al. 1999) at A 2000 A in the rest frame 
(assuming Ho = 65km s~^Mpc~^, an open cosmology with Q = 0.3, and zero extinction). Broken curve: eq. (1) 
with Lo — 2 X 10*^ ergs s~^ and to — 20 s. This underestimates the energy in the optical transient of GRB 990123 
optical transient. Solid curve: eq. (1) with Lo = 2.5 x 10*** ergs s^^ and to — 10 s, which we adopt for the optical 
transient of a typical GRB. 



2. The Optical-UV Flash 

Let L = AirdP/dVl where dP/dVt is the radiated power per solid angle in our direction. It 
is not yet known what should be taken for the "typical" L{t) for a GRB at optical - soft X-ray 
wavelengths. 

In Figure 1 we show vLy{t) of GRB 990123, at redshift z = 1.6,^ at rest-frame wavelength 
~ 2000 A. 



Assuming Ho = 65 km s ^ Mpc ^ , an open cosmology with Q, = 0.3, and zero extinction. 
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As a simple fitting function to describe the time-dependent luminosity, we adopt 



'[l + (Vio)2]2 Vl3.6eV 

A spectral index (3 = —0.5 is suggested by simple models for the reverse shock (Waxman &; Draine 
2000). Steeper spectral indices have been observed for optical-UV afterglows: for GRB 010222, at 
redshift z > 1.475 (Masetti et al. 2001) Lee et al. (2001) find P = -0.90 at (1 + z)t = 1.7 x 10^ s, 
Stanek et al. (2001) find /? = -1.07 at {l+z)t = 1.7x 10^ s, and Masetti et al. (2001) find /? = -0.55 
for 8 X 10^ < (1 + z)t < 1.6 X 10^ s. 

With Lq = 2 X 10^^ ergs s~^, to = 20s and /? —0.5, eq. (1) approximates the light curve 
of GRB 990123 at 50 < t < 500 s, although falling at least a factor 5 below the observed peak 
bightness at t = 20 s (see Fig. 1). 

The Ljy oc dependence of equation (1) at t ^ is required to approximate the observed 
rapid decline in brightness of GRB 990123 between 20 and 300 sec (sec Fig. 1). However, at later 
times afterglows appear to fade less rapidly: for GRB 010222, oc t-0-75±0.05 u ^ iqA < 
{l + z)t<3x s, and oc t'^-^ for 8 x 10"^ < (1 + z)t < 3 x 10^ s (Stanek et al. 2001; Masetti 
et al. 2001). For the present study we are concerned with the lightcurve at t < 10^ s, which we will 
approximate using equation (1). 

GRB 990123 was in the top 0.4% in gamma-ray fiuence for GRBs detected by the BATSE 
instrument, and LOTIS observations further suggest that most GRBs detected by BATSE have 
lower ratios of optical to gamma-ray emission (Williams et al. 1999). Accordingly, we adopt Lq = 
2.5 x 10^^ ergss"^ - i.e., a peak luminosity about 1/40 of what was observed for GRB 990123 - and 
to ~ 10 s as "typical" parameters. For /? = —0.5 this corresponds to the optical transient having an 
energy between 1 and 13.6 eV given by E^oT = 27rLoto[l- (l/13.6)(l-|-/3)]/(l+/3) = 1.57 x 10^° ergs; 
for ^ = -1, EoT = 27rLoto In 13.6 = 4.1 x 10^° ergs. 



3. Radiative Transfer 

We divide the cloud into spherical shells j = 1,...,N oi uniform thickness AR, with outer radii 
Rj and midpoints Rj = Rj — 0.5Ai?. At each radius r, let tr = t — r/c be the "retarded" time. 

The shell thickness AR is small enough [njiAR < 10^^ cm~^) so that an individual shell 
is optically thin in dust, but the shells may be optically thick to photoionizing radiation, or in 
absorption lines of H2. We therefore calculate rates for photoionization, photodissociation, or 
photoabsorption averaged over the shell volume 



AVj = (47r/3) [R'^ - {Rj - ARf] 



(2) 
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At frequency v, the optical depth contributed by shell j is 

ATjiiy) = dTdjW) + <5th2j(z^) + + + 5t^^oj{v) + 8T^^+^j{v) + ^ Jr^j (3) 

a 

where ^t^j, Stu^j, j-: ^"^nj^ ^'^ne°,j^ "^"^He+j contributions in shell j due to ex- 

tinction by dust, and continuous absorption by H2, Hg", H, He°, and He"*", respectively; Taj is the 
contribution of an individual H2 absorption line a (see §3.3). Absorption by metals not in dust 
grains (initially O, N, Ne, and S will be most important, and C, Mg, Si, and Fe will enter the 
gas phase in regions where grains are vaporized) has been neglected in eq. (3), because the present 
calculation is limited to the gas phase species H, H"*", H2, H^, He, He"'", and He^"*". We have thereby 
underestimated the attenuation at hu > lOOeV. 

We divide the electromagnetic spectrum into 3 intervals: hv < 11.1 eV, 11.1 < hu < 13.6 eV, 
and 13.6 < hu < 20keV. For purposes of modelling the heating of dust grains by the UV-optical 
radiation, we include the contribution of H2 absorption lines only within the interval 11.1 < hv < 
13.6 eV; the summation in equation (3) is restricted to lines a within this interval, with (5tq 
calculated from eq. (14) below. 

The total optical depth from r = to Rj is 

3 

Tj(y,tr) = ^ATi{v,tr) . (4) 

i=l 



3.1. Absorption and Scattering by Dust 

The contribution by dust in shell j is 

(^Tdj(i/) =ndj7ra^(Qabs(z^) + Qsca(z^))A-R , (5) 

where aj{tr) is the dust radius in shell j at retarded time tr, and (5abs(i^) and Qscal'^) are the usual 
efficiency factors for absorption and scattering. The radiation from the GRB is concentrated in a 
shell-like "pulse" with a spatial thickness of ~ 10 light-sec (see Fig. 1). Scattering is important 
because it reduces the intensity in the pulse; the scattered photons will fall behind the outward- 
propagating pulse unless the scattering angle is very small. For scattering angle ^ <C 1, after 
travelling a distance 6R the scattered photon will fall behind the unscattered pulse by a time 

sne^ , , SR f e Y 

Stf^— — = 1.4s— -] , (6) 

2c 10^® cm \arcmmy 

so that the time delay will be significant for scattering angle 6 > Sarcmin. The dust grains will have 
radii a ~ 3 x 10~^ cm (see §5.1). For a = 0.25 fim, the characteristic scattering angle is ~ 3 arcmin 
at 2 keV, and ~ 10 arcmin at 0.5 keV (see Figs. 3 and 4 of Smith & Dwek 1998) so we take the 
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scattering efficiency factor to be 





Qsca(i^) = S 1 





hi^ <leV , 

leY <hu <1 keV 

hu > 1 keV . 



In Fig. 2 we show Qabs calculated for a = 0.3;^m graphite and silicate grains. At hi^ > 1 keV, the 
silicate absorption efficiency is much larger because of the inner shell absorptions from Mg, Si, and 
Fe. In the present calculation we do not distinguish between grain types, and we need a simple 
estimate for the absorption efficiency which can be rapidly evaluated. As a simple approximation, 
we adopt 



which can be seen in Fig. 2 to approximate the average absorption cross sections of graphite and 
silicate grains. 



hi/(eV) 

Fig. 2. — Absorption efficiency Qabs for a — O.Sfim radius grapfiite and silicate grains. 

When grain destruction has reduced the grain radius to a ^ 3 x 10^^ cm, eq. (8) will overesti- 
mate Qabsi when this is true, however, grain destruction is nearly complete and the error in overall 
grain destruction will be unimportant. 




T-rry 




eq. (8) 
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3.2. Continuous Absorption by Atoms, Molecules, and Ions 

Continuous absorption by species X contributes 

5Tx,j(y) = n{Xj)ARax{y) , (9) 

where n{Xj,tr) is the abundance of X in shell j at retarded time tr, and crx{i^) is the continuous 
absorption cross section. 

The probability that a photon which enters the shell will be absorbed or scattered is (1— e"^'^^ ). 
For a homogeneous shell, a fraction Atx^/^tj = nx,jcrx^R/^Tj of the absorptions will be due 
to continuous absorption by species X (e.g., X + hv ^ X+ + e~). The photon absorption rate per 
X, averaged over the shell, is 

Integration over frequency v is accomplished numerically, with Atj and Tj-i evaluated using equa- 
tions (3) and (4). 



3.3. Line Absorption 

Line absorption is important for H2. Consider a specific absorption line a, with central fre- 
quency We assume that the H2 has a MaxwcUian velocity distribution with broadening param- 
eter 6, and we compute the dimensionless equivalent width 

/I 
^ {1 - exp [-N^{Rj)c7o,{v)]} , (11) 

where Uai^') is the (Voigt profile) absorption cross section for absorption line a, and 

3 

NaiRj) = I] {^2{Va, Ja))]i AR (12) 
1=1 

is the column density of the vibration-rotation state (?;„, Jq,) producing absorption line a. 
The incremental equivalent width contributed by shell j is 

SW^j = WaiRj) - W^{Rj-i) . (13) 

If we approximate the line absorption as due to many lines randomly distributed over the frequency 
interval [t'min, z^max], the optical depth contribution due to line a averaged over the wavelength 
interval is (see Draine &; Bertoldi 1996) 

= 1 r T for Va e [t'min, I'max] (14) 

Inifmax/l'min] 
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and the mean photoabsorption rate in radial shell j for excitation out of level Xa by the absorption 
line a is 

47rR]N^{Rj)h ■ ^'^^ 

Strong line absorption by H2 occurs primarily in [/ifmin, ^z^max] = [11-2, 13.6] eV, so we include the 
optical depth correction (14) only for this interval. 



4. Physical Processes 

4.1. Dust Destruction by Thermal Sublimation 

The thermal sublimation rate from a grain at temperature T can be approximated by (Guhathakurta 
& Draine 1989, hereafter GD89) 

1 

da I 



where m is the mean atomic mass, and B is the chemical binding energy per atom. We take 
I'o = 1 X lO^^s^^, B/k = 7 X lO^K, and p/m = 1 x lO^'^^cm"^ as representative values for refractory 
grains (GD89; WDOO). We assume that all of the energy of absorbed photons (including X-rays) 
is converted to heat - we neglect the energy of photoelectrons and X-ray fluorescence, but this is a 
smal fraction of the total absorbed energy. 

In the middle of shell j, the grain temperature T is determined by (WDOO) 

r exp (-r,j_i - 0.5At,j) QabsT^a'' = {Q)T^TTa''aT^ - Aira^^-^B , (17) 

where {Q)t is the Plank-averaged absorption efficiency. For the temperature range of interest for 
dust sublimation, 2000 K <T< 3000K, we approximate (WDOO) 

^ 0.1(a/10-5em)(r/2300K) 
^^^^ ^ l + 0.1(a/10-5cm)(T/2300K) ' ^ ' 

intermediate between the emissivities of astronomical silicate and graphite. 



4.2. Dust Destruction by Electrostatic Stresses? 

Waxman & Draine (2000) noted that a highly-charged dust grain can also be destroyed by 
"Coulomb explosion" if the electrostatic stress exceeds the tensile strength of the grain materials, 
when the grain potential reaches 

U{,^ ^ 3000V(Sm£«/10^^ dyncm-2) 1/2 (a/10-5 cm.). (19) 
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The importance of this process is controversial. Waxman &; Draine argued that submicron grains 
would likely have high enough tensile strengths (S'max > 1 x lO^^dyncm"^) that highly-charged 
grains would not undergo fission by this process, but would instead be gradually eroded by the 
process of "ion field emission", where singly-charged ions are emitted one-by-one when the grain 
potential exceeds a critical value 

UiFE ~ 3 X lO^V (a/10-^ cm) (20) 

(MuUer & Tsong 1969; Draine & Salpeter 1979). Pruchter, Krolik, & Rhoads (2001), on the other 
hand, have argued for lower tensile strengths (S'max ^ lO^dyncm"^), in which case grain fission 
would take place when the grain potential reached C/frac ~ 300V(a/10~^ cm)(5max/10^ dyncm"^)^/^. 
Fruchter et al argue that "the large flux of energetic photons bombarding the grain is likely to dam- 
age the grain's crystalline structure" , reducing the tensile strength. 

While chemical bonds will undoubtedly be disrupted by ionization, it seems likely that chemical 
bonds will be promptly reestablished in the warm grain - i.e., "annealing" will take place. In this 
case, a grain which has already been charged to the ion-field emission threshold potential ?7ife 
would emit one ion per escaping photoelectron. To sustain ion field emission, the grain must 
emit photoelectrons with sufficient energy to overcome the potential [/ife = 3000y(a/10~^ cm). 
Consider, for example, a grain with radius a = 3 x 10~^ cm, charged to C/jfe ~ 9 kV: only 
photoelectrons with kinetic energy at the grain surface E > eU ~ 9keV are able to escape. For 
a grain with composition MgFeSi04, a photon of energy hf = lOkeV (for example) has a total 
absorption coefficient a = 600 cm~^ in the grain material, but only ~ 4% of the absorptions produce 
photoelectrons with E > 9keV,^ for an effective absorption coefficient 24cm~^. The total electron 
density in the grain is 1 x 10^^ cm~^, so the effective photoabsorption cross section per electron is 
~ 24 X 10~^^ cm~^ - a factor of 24 larger than estimated by Waxman & Draine (2000), but a factor 
~ 10'^ smaller than the absorption cross section just above the Si K edge which Fruchter et al used 
to estimate the rate of photoelectric emission driving grain fragmentation by electrostatic stresses. 

Grain destruction by ion field emission requires a fluence -Fife ~ {p/i^)/ctiFE, where p/m ~ 
10^^ cm""^ is the atomic density and oife is the effective absorption coefficient for producing photo- 
electrons with energies > cUife- For C/jfe = 9kV, we have aiFE = 24cm~-^ and Fife = 4x 10^^ cm~^ 
of hu > lOkeV photons. For /? = —0.5 and distance D, the fluence of hu > lOkeV photons is 
2.2 x 10^-'^(Loto/10^^'^ergs)(pc/F))^cm~^ (neglecting absorption), so we might expect grain de- 
struction by ion field emission at distances D < Ipc. However, as we will see below (§6), grain 
destruction by thermal sublimation is effective out to distances of several pc. 



^hu = lOkeV photons absorbed by the Mg, Si, and Fe K shells (with photoelectric thresholds Et = 1.32, 1.87, 
and 7.11 keV) produce photoelectrons with hv — Et = 8.68, 8.13 and 2.89 keV, respectively. 
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4.3. Photoexcitation of H2 

Let index i = 1, ...,N denote the vibration-rotation levels of the electronic ground state X^5]+. 
At a given point in the cloud, let pi(t) = 2n(H.2{vi, Jj))/nH be the fraction of H nuclei which are in 
vibration-rotation level i. 

Prior to the GRB, H2 is almost entirely in the first two or three rotational levels (J=0,l,2) 
of the v=0 vibrational level of the electronic ground state (X^S+). The first UV photons to 
be absorbed will photoexcite the H2 to the B^S+ and C^n^ electronic states (via Lyman or 
Werner band transitions), which requires hu > 11.2 and 12.3 eV, respectively) or photoionize the 
H2 to H^, which requires hv > 15.4 cV [for H2(v = 0)]. H2 which is photoexcited to the B or C 
states will decay back to the ground electronic state in ~ 10~^ s, but typically to a vibrationally- 
excited level {e.g.,v = 5) with the rotational quantum number changed by A J = 0, ±2. The 
lifetimes of the vibrationally-excited levels are long compared to the timescale for photoexcitation 
or photoionization, so depopulation of the vibrationally-excited levels of H2 will be primarily by 
UV photoexcitation and photoionization. 

The intensity of the radiation field at lOOOA relative to the local interstellar radiation field 
estimate of Habing (1968) is 



A. 



4 X 10~^^ ergs cm""^ 10^^ ergs s 

Stimulated emission in the UV transitions is negligible if x ^ lO"*^^. We therefore assume 
that photoexcitation out of level i of to vibration-rotation states of B^S+ and C^II^ will 

be followed (immediately) by spontaneous decay either to bound levels j of the ground electronic 
state X^S+, or else to the vibrational continuum of X^S+ (i.e., photodissociation) . H2 formation 
on grains or in the gas will be neglected since we will be concerned with conditions where the 
H2 is destroyed in a matter of seconds or minutes. Collisional deexcitation can be neglected if 
nn/cm^^ <^ X- The spontaneous decay via quadrupole transitions out of vibrationally-excited 
level i can also be neglected if x ^ 10^- 

Let Tji be the rate for photopumping of H2 out of level i into level j, and let cf^*^^ ^"^^ C!f^^ be 
the rates for photodissociation and photoionization out of level i. If 



-'-II — 



(22) 



then the H2 abundances pi = 2n[H2(fi, Ji)]/nH evolve according to 



^=j:Tijpj- (23) 

j 

Following Draine &; Bertoldi (1996), we consider the 299 bound states of H2 with J < 29, and 
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calculate equivalent widths Wa using Lyman and Werner band oscillator strengths from Abgrall et 
al. (1993a,b) and Roueff (1993). 

Let a denote a Lyman or Werner band transition from a specific vibration rotation state i 

of the electronic ground state to an excited state u{a). Photocxcitation rates (a in shell j are 
calculated using eq. (15). Let Pj^u be the probability that clcctronically-excitcd state u will decay 
to vibration-rotation state j of the electronic ground state. Averaged over the shell, the rate for 
UV pumping from vibration-rotation level i to vibration-rotation level j is 

Tji = Y^Pj,u{a)<a (24) 

a 

where the summation is restricted to transitions a out of level i. The photodissociation rate out of 
level i is 

=Y^mssHa)\ (25) 

a 

where paiss ("") is the probability that level u will decay by a transition to the vibrational continuum. 
We neglect line absorption by the ion, as it is likely to be of secondary importance. 



4.4. Yi2 + hv + e- 

The different vibrationally-rotation levels of the electronic ground state of H2 will have different 
energy thresholds for photoionization. We will, however, neglect these differences and assume that 
all H2 in the ground electronic state requires hv > 15.4 eV to be photoionized, with the cross section 
as given by eq.(17-19) of Yan et al. (1998). The photoionization rate (^^^ averaged over the shell 
is evaluated using eq. (10). 



4.5. B.+ + hu + R 



has 19 vibrational levels; for photoionization of H2(f = 0) the vibrational distribution 
is expected to peak at v = 2, with 80% of the population in v < 5 (Busch Sz Dunn 1972). With 
radiative lifetimes ~ 10^ s, there will be negligible deexcitation by spontaneous decay during the 
optical transient. Averaged over the vibrational distribution, the photodissociation cross section 
found by von Busch &; Dunn (1972) can be fit by 

<, = 2.7xl0--cm^(^)'(l-^)' torft.<29eV. (26) 

giving a peak cross section cr 3 x 10~^^ cm^ at hu 7.3 eV (A ^ 1700 A). The photodissociation 
rate (^^^ averaged over the shell is evaluated using eq. (10). 
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4.6. Rt + hiy^ 2H+ + e" 



For photoionization of we adopt a cross section 



(7 = 9.3 X 10"^^ cm^ 



hi' 



15.4 eV 



for hv > 15.4 eV 



(27) 



which approximates the photoionization cross section calculated by Bates & Opik (1968) for the 
ground vibrational state of H2'(u = 0). The photoionization rate (^^^ averaged over the shell is 
evaluated using eq. (10). 



4.7. H + /ii/ ^ H+ + , 



Atomic H from the photodissociation of H2 or H2' can be photoionized by photons with hu > /h 
with a cross section (Osterbrock 1989) 



a = 6.30 X 10~^^ cm^ 



1.34 



The photoionization rate Ch**^ averaged over the shell is evaluated using eq. (10). 



hu 



13.6 eV 



-2.99 



0.34 



hu 



13.6 eV 



-3.99 



for /iz/ > 13.6 eV. (28) 



4.8. Re + hu ^ He+ + e" and He+ + hu ^ He^+ + e" 
The cross sections for photoionization of He and He+ are (Osterbrock 1989) 



a = 7.83 X 10"^^ cm^ 



a = 1.58 X IQ-^^cm^ 



1.66 



1.34 



hu 



24.6 eV 



hu 



-2.05 



-2.99 



0.66 



0.34 



hu 



24.6 eV 



hu 



-3.05 



-3.99' 



for /iz/ > 24.6 eV (29) 



for hu > 54.4 eV (30) 



54.4 eVy V54.4eV^ 
The photoionization rates <^^^' and C^jj averaged over the shell are evaluated using eq. (10). 



5. Photochemical Evolution 



5.1. Initial Conditions 



As initial conditions we take He to be neutral, and 100% of the hydrogen in H2, divided equally 
between the first two rotation levels: 



^ n[R2{v = 0,J = 0)] _^ n[il2{v = 0,J = l)] 



(31) 
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rig. o. — Snapshots of the ionization/dissociation/dust destruction fronts at three values of the retarded time tr 
for our standard GRB [L^ oc given by eq. (1) with Lq = 2.5 x 10^*ergss~^, to = 10s, P = —0.5] in a cloud 

with riH = 10"^ cm""^. Open triangles: dust mass (a/ai)^; circles: 2n(H2)/nH; open diamonds: 10 x n{¥L)/nu; filled 
diamonds: n(H+)/nH. 



We assume that dust grains of initial radius Oj = 3 x 10 ^ cm initially contribute a mass equal 
to 1% of the H mass, with initial number density 



rid 



O-Olnnm-H 

(4vr/3)pa3 



(32) 



where p ~ 2 g cm ^ is the assumed density of the grain material, nn is the number density of H 
nuclei, and mn is the mass of one H atom. Thus the initial geometric cross section per H atom is 



^7ra2 = 2.1 X 10-22 cm^H-i 



(33) 



Since the radiation propagates only in the outward direction (we are neglecting scattered 
light), the time evolution of each shell depends only on interior shells. We therefore calculate the 
complete time evolution for each shell, beginning with shell 1 and proceeding outward. For each 
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shell j, we begin integrating at retarded time U = 0. For each value of U we use stored information 
for interior shells i < j to compute the time- and frequency-dependent optical depth Tj_i due to 
material interior to shell j. 

The calculation for shell j is terminated when tr > lOOio = 10^ s [see eq. (1)] and nearly all 
of the flash energy has been expended, or if all of the following conditions are fulfilled: 1. all H2 
has been photodissociated to 2H and then photoionized into 2H"'", or photoionized first to and 
then to 2H+ [our criteria for this are J^iPii^i^ Ji) < 10~^> ^(H) < 10"^ and x(H+) < 10"^]; 2. the 
dust has been destroyed (a = 0); 3. helium is fully photoionized to He^"*". 

5.2. Integration Scheme 

A very simple explicit scheme is used to advance the abundances in each shell: 

y{t + At)=y{t) + {dy/dt)tAt , (34) 
where y is a, n(H), n(H+), n(Re^), n(Re'^), n(He^^), or H2 vibration-rotation level population pi. 
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For all cases we used radial zone thickness AR = 6 x 10^^ cm. At a given distance r and 
retarded time U, we take the time step At to be 



t^^ = max 



At = mm[eatd, f-ptp] 

_rn /-(pi) //-ipi) j_ /-(j>d)\ Api) /-(pi) 
-'mAh ' v'»h+ '''^He ''»He+' 



^-1 



(pi) n(H2) 
n(H+) 



(35) 
(36) 

(37) 



td, is the shortest time scale for destruction, while tp is the timescale for production of H^. The 
rates Tu arc given by equation (22). Wc used {ed,ep) = (0.2,0.2) for all cases; we verified that the 
computed results were insensitive to factor-of-two variations in and e^. 



6. Results 

Radial abundance profiles are shown in Fig. 3, for three values of the retarded time tr, for an 
optical transient with the light curve given by eq. (1) with spectral index /3 = —0.5, located in a 
molecular cloud of density rifj = 10'^cm~^. 

In Fig. 4 wc show abundance profiles for the same three values of retarded time, for the same 
optical transient but situated in a cloud of density nn = lO^cm"^. In its initial molecular form, 
an individual shell is opaque to ionizing radiation hv < 70 eV, and is highly opaque in individual 
H2 resononance lines; it is therefore essential to use the shell-averaged rates given by eq. (10) for 
continuum absorption, or (15) for line absorption. As one would expect, the increased opacity of 
the medium delays propagation of the destruction fronts. 

In Figs. 5 and 6 wc show the locations of the dust destruction, dissociation, and ionization 
fronts as functions of time for our standard optical transient in clouds of density njj = 10^ and 
lO^cm"^, and we also show the front locations as functions of time for a GRB with a softer UV 
spectrum characterized by /3 = — 1. 



6.1. Dust Destruction 

In Fig. 3 we see that by = 2s the dust destruction front has already reached R 3.7 x 
10^^ cm, significantly ahead of the ionization/dissociation front at ~ 1.9 x 10^^ cm. At late times 
the dust destruction front "stalls" when the transient luminosity has dropped to the point where 
it can no longer heat dust at this distance to the T « 3000 K temperatures required for significant 
sublimation, while the ionization/dissociation front continues to advance at a rate proportional to 
the luminosity. As a result, for this case the ionization/dissociation front catches up to the dust 
destruction front (at i « 70 s) and moves past it at late times (see Fig. 5). The dust destruction front 
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Fig. 5. — Location of the dust destruction front, H2 destruction front, and ionization front as functions of time for 
our model optical transient in molecular gas with nn = 10^ cm~'^, for /3 = —0.5 and -1. The ionization front and H2 
destruction front are nearly coincident. 



(see Fig. 3) has a thickness 5R ~ 1 x 10^^ cm, corresponding to a gas column density ~ 10^^ cm ^, 
and an optical-uv extinction (prior to dust destruction) 0.4. 

In Figs. 5 and 6 we show the location of the dust destruction front as a function of retarded 
time tr for our adopted optical transient (for spectral indices /? = — 1 and /? = —0.5) and cloud 
densities nn = 10^ and lO^cm"^. For both cases we confirm the conclusion of WDOO: the dust 
destruction front initially advances well ahead of the ionization/dissociation front, with the dust 
heated primarily by hv < 11 eV radiation which is not absorbed by H, He, or H2(v = 0). For 
nil = 10^ cm^'^ and (3 = —0.5 the dust destruction radius ~ 1.2 x 10^^ cm for our adopted 
optical transient, while for lO^cm"^ we find ~ 7.7 x 10^^ cm. 

For a transient with peak luminosity 3 x 10'^^ ergs s~^ in l-7eV photons - close to the peak 
1-7.5 eV luminosity 2.4 x 10^^ ergs s^^ for our adopted light curve with (3 = —0.5 - WDOO estimated 
the grain destruction radius to be Rd ~ 1.4 x 10^^ cm for nn < 10^ cm~^, and 1.3 x 10^^ cm for 
nu = lO'^cm"'^. The grain destruction radius estimated by WDOO for these two cases is in good 
agreement with our detailed modelling. 
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Fig. 6. — As in Fig. 5, but for cloud density nn — 10'' cm~^. 

6.2. Photodissociation and Photoionization 

In Fig. 3 we can see that the dissociation front and ionization front are merged (i.e., they 
do not separate). However, the ionization front is relatively sharp (with the fractional ionization 
varying from 0.1 to 0.9 over a column density AA'h ~ lO^^cm"^), whereas the dissociation front 
(zone where H atoms are present) tends to be broader, and to extend ahead of the ionization front. 
For nil = 10^ cm~^ (see Fig. 3) the dissociation front extends ~ 1 x 10^'^ cm, corresponding to the 
column density nuSR ~ 10^'^ cm~^ required for H2 to absorb most of the photons between 11.1 
and 13.6 eV (Draine & Bertoldi 1996). As expected, the ionization front is unresolved, since the 
column density of a single shell (5 x 10^^ cm~^) corresponds to an optical depth of ~ 35 at the 
Lyman edge [hence the importance of using eq. (10) to obtain the average rate of photoionization 
in an optically-thick shell]. This is even more pronounced for the nn = lO^cm"^ case, for which 
we have used the same shell thickness AR as for nn = 10^ cm~'^. 



6.3. Vibrationally-Excited H2 

Following absorption of a Lyman or Werner band photon, about 85% of the time an H2 molecule 
will undergo a spontaneous radiative decay (in ~ 10~^ s) to a vibrationally-excited level of the 
ground electronic state. The vibrationally-excited H2 will persist until it is either photoionized 
or absorbs another Lyman or Werner band photon (on the ~ 10^ s timescales of interest here. 
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Fig. 7. — Column densities of and vibrationally-excited H2 for tih = 10^ cm~^ and 10* cm~^, and [3 = —1 and 
-0.5. The oscillations in A*'(Hj) are due to inaccuracies in the numerical method as the ionization front traverses 
a radial zone. The column density of quickly saturates at ~ 3 — 10 x lO^^cm"^, but the column density of 
vibrationally-excited H2 continues to rise until the ultraviolet flash fades. 

spontaneous emission in the quadrupole lines is negligible, as are collisional processes). In Figure 7 
we show the column density A^(H2) of vibrationally-excited H2 as a function of time. For all cases 
considered, A^(H2) increases until the UV irradiation ceases [unlike the H^, which quickly stabilizes 
at a value A^(H^) ~ 10^^ '' cm~^]. Photoionization acts to reduce the amount of vibrationally- 
excited H2 present at any instant, so GRB with softer UV-EUV spectra would be expected to have 
larger column densities of vibrationally-excited H2. Our calculations confirm this: Fig. 7 shows that 
an optical transient with spectral index /? = —1 (i.e., L^, oc produces more vibrationally-excited 
H2 than a burst with /? = —0.5. 

As we see below, the large column densities of vibrationally-excited H2 will produce strong 
absorption lines. 

6.4. Emergent Spectrum 

The radiation reaching us from the optical transient is filtered through the dust and gas. Unless 
and until the photoionization front reaches the edge of the gas cloud, photoelectric absorption 
by H will impose a cutoff at 13.6 eV, with photoelectric absorption by H2 further contributing 
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Fig. 8. — Dust-attenuated light curves at rest wavelength Ao = 1800 A for different values of the cloud radius -Rmax, 
for a GRB with (3 = —0.5 in a cloud with tih = 10^ cm~"^. A dust extinction cross section/H of 4.2 x 10~^^ cm^ is 
assumed. For a GRB a,i z — 1, the U band apparent magnitude scale is indicated at right. 

to absorption above 15.4 eV. At X-ray energies hv > 10 keV the gas cloud may again become 
transparent, allowing the X-ray flash and afterglow to be observed, but our attention here will be 
limited to energies hi' < 13.6 eV, where there will be four main contributions to absorption: dust, 
the Lyman and Werner bands of H2, photodissociation of H^, and the Lyman lines of H. 

The extinction by dust will aff'ect whether the optical transient is observable. In principle there 
could be spectral features (e.g., the 2175 A feature) in this absorption, but our understanding of the 
optical properties of dust in these dense regions is such that no realistic estimate is possible. In Fig. 
8 we show the emergent luminosity after extinction by dust for ambient density nn = lO^cm"^. 
Light curves are shown for various assumed cloud radii; for each case the initial visual extinction Ay 
to the GRB is indicated. In the case of i? = 1 x 10^^ cm cloud radius, the ionization/dissociation 
front reaches the cloud surface at t ~ 30 s; after this time, there will of course no longer be 
absorption features of H2 appearing in the spectrum of the GRB. For R > 1.4 x 10^^ cm, however, 
the ionization/dissociation front is contained within the cloud, and the II2 absorption lines will 
persist in the spectrum of the GRB afterglow at late times. 

For purposes of discussion we focus on our standard lightcurve with /? = —0.5 in a cloud of 
density nn = lO^cm"^. By destroying the dust out to ~ 1.2 x 10^^ cm, the GRB has reduced 
the extinction between us and the GRB by AA 5 mag. Thus if, for example, the dust in the 
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Fig. 9. — Transmission spectrum of the intervening H2, H, and at selected times, for a GRB in a cloud with 
nu ~ 10^ cm~^, as observed by a spectrograph with R — A/FWHMa = 350. The GRB is assumed to have a spectrum 
with P = —0.5 (i.e., L^, oc i/""'^). Continuum absorption by dust contributes additional attenuation. At resolution 
R = 350 a number of A > 1110 A absorption "features" from vibrationally-excited H2 are evident, and are labelled 
by wavelength. Some of the stronger lines contributing to these features are listed in Table 1. 

cloud had initially presented an extinction of, say, 8 mag to the GRB, the GRB afterglow would 
be subject to only 3 mag of foreground dust. 

The column densities tend to be considerably smaller than those of H2 (see Figure 7), so we 
have not attempted to calculate the column densities of the individual vibration-rotation levels and 
the —5- C^Hu ultraviolet absorption lines out of these levels. Most of the electronic absorp- 

tion by will be in X'^'^^ — ^'^'^t photodissociating transitions, which contribute continuous 
absorption which we have approximated by eq. (26). 

In Fig. 9 we show the transmission of the foreground medium for the nn = 10^ cm~^ cloud and 
spectral index f3 = —0.5, for spectral resolution R = 350 characteristic of the grism on the Swift 
UVOT instrument. Fig. 9 shows only the absorption due to lines of H2 and H, and the continuous 
absorption of H^. Voigt profiles were used for both the H and H2 lines. Over the < 10^ s timescales 
of interest here, the H and H2 gas is essentially collisionless. For the H2 lines we have therefore used 
a Doppler broadening parameter b = 3 km s~^ characteristic of "microturbulence" in the quiescent 
pre-GRB molecular gas. The H atoms, on the other hand, will have a velocity distribution resulting 
from photodissociation of H2 and H^. The mean kinetic energy per H atom is ~ 0.15 eV (Stephens 
& Dalgarno 1973) corresponding to a Maxwellian distribution with b ~ 4 km s^^; this will add in 
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The A > 1110 A spectrum for the same conditions as Fig. 9 at spectral resolutions R = 100, 350, and 



quadrature to the microturbulence, so we take b = 5 km s for the H atoms. 

The spectrum consists of hterahy thousands of saturated narrow hues, but even at i? = 350 
resolution, a number of absorption features are conspicuous. In Figure 9 we have labelled 32 
conspicuous features longward of A = 1111 A. These features are ah blends of a number of nearby 
lines. For each feature in Table 1 we list the three strongest lines within ±X/{2R) = ±1.3 x 10~^A of 
the central wavelength of each feature. The strongest lines tend to be saturated (e.g., the strongest 
line in the conspicuous 1277. 3A feature is Lyman 8-6R(l) 1277. 39A with a central optical depth 
To = 1150 and an equivalent width Wx/X = 8.3 x 10^"'). 

Figure 10 shows the transmission spectrum observed with three different spectral resolutions: 
R = 100, 350, and 1000. While the R = 1000 spectrum of course shows the most structure, we see 
that the R = 350 resolution is quite well- matched to the widths of the prominent absorption blends. 
If the resolution is degraded to i? = 100, most of the absorption blends still appear as well-defined 
minima, though a few do not (e.g., A = 1621. 7A.) For R = 100 observations, it would appear that 
a signal-to-noise ratio of 10 or more would be sufficient to test for the presence or absence of the 
stronger absorption features (e.g., 1607.0, 1460.0, 1401.9, 1277.3, 1254.8 A). 

Because the lines are saturated, the overall spectrum is not sensitive to the precise amount of 
vibrationally-excited II2. To see this, in Figure 11 we show the R = 350 spectrum for three different 
values of the column density A^(Il2) of vibrationally-excited II2, ranging from A^(H2) = 10^®'^^ cm~^ 
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(the value at t = 2s for the f3 = -0.5, nu = 10^ cm-^ case) to iV(H^) = lO^^-^" cm-^ (the value 
at t = 10^ s for the /3 = —0.5, tt-h = 10^ cm~^ case). While the depth of the absorption features 
does increase with increasing A^(H2), the increase is fairly modest despite the order-of-magnitude 
change in A^(H2). 



7. Observational Prospects: HETE and Svi^ift 



The presence of absence of the absorption features in Table 1 will indicate whether or not 
GRBs occur in the vicinity of molecular gas. The present paper has carried out simulations of the 
effect which the optical and ultraviolet radiation from a GRB will have on nearby gas and dust, 
assuming that the GRB progenitor is located in a uniform molecular medium of density nn = 10^ or 
10"^ cm~^. It is important to realize, however, that vibrationally-excited H2 can be produced even if 
the GRB is in a low-density region, but there is molecular gas on the line-of-sight within a distance 
of a few parsecs from the GRB. Absorption features will be conspicuous for vibrationally-excited 
H2 column densities A^(H2) > 10^^ cm^^. This will be produced provided the total H2 column 
exceeds 10^^ cm~^, and the GRB flash includes a flux (l/47ri?^) J dt[i'L^/{hi>)]^^^^^ > Fmin, where 
-^min ~ 5 X 10 cm" of photons at A ?a 1000 A. For the light curve (1), this flux occurs at a 
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Table 1: A > lllOA Absorption Features from UV-pumped H2 



A(A)« 



Three strongest contributing lines 



1704.7 
1669.3 
1654.6 
1637.5 
1621.7 
1607.0 
1590.5 
1578.6 
1545.8 
1517.5 
1504.2 
1494.6 
1487.1 
1460.0 
1436.2 
1401.9 
1394.8 
1374.8 
1357.5 
1344.6 
1334.5 
1324.3 
1306.0 
1277.3 
1254.8 
1199.2 
1180.6 
1162.3 
1145.6 
1134.6 
1128.4 
lllG.O 



3- 13P(l) 1704.51 
2-llR(l) 1667.46 

1- lOR(l) 1654.99 

4- 12R(l) 1638.90 

2- 10R(l) 1620.74 

6- 13P(l) 1607.50 

3- 10R(l) 1588.76 

7- 13R(l) 1577.02 



1-8R(1 

0- 7R(l 
10-13P 

1- 7P(2 

1- 7R(1 
0-6P(2 

3- 7P(l 

2- 6R(l 
0-5R(l 
10-9R( 

4- 6R(l 

2- 5R(l 
0-4R(l 

3- 5R(l 
9-7R(l 

0- 3P(l 

1- 3R(1 

1- 2R(1 

2- 2R(l 
1-4Q(1 

1- 1R(1 

2- lP(3 
2-lR(l 
1 3Q(1 



1544.90 
1516.23 
1) 1504.15 
1491.71 
1486.63 
1460.17 
1435.05 
1402.14 
1393.97 
) 1374.49 
1355.56 
1345.40 
1333.80 
1323.28 
1306.21 
1276.82 
1253.94 
1198.49 
1180.42 
" 1161.29 
1144.71 
1135.36 
1128.22 
lllG.ol 



3- 13R(0) 1702.63 

4- 13R(l) 1667.28 

5- 14P(l) 1654.19 
3 llP(l) 1636.33 

2- 10P(l) 1623.51 
5-12R(0) 1608.43 

3- 10P(l) 1591.31 



0- 8R(l 

1- 8P(1 

2- 8P(l 
10-13P 
1-7P(3 
1-7P(1 

4- 8R(l 

5- 8R(l 
0-5P(3 

0- 5P(l 

1- 5P(2 
4-6R(2 

2- 5R(0 
0-4R(0 

3- 5R(0 
3-9R(l 

0- 3R(2 

1- 3R(0 

1- 2R(0 

2- 2P(l 

1R(1 
0-3Q(l 
2-4Q(2 
2-lR(0 

1 3R(1 



1576.88 
1547.54 
1517.44 

2) 1504.92 
1495.22 
1489.08 
1460.77 
1436.09 
1402.66 
1396.23 
1373.66 
1357.35 
1345.04 
1333.48 
1322.80 
1305.59 
1276.33 
1253.50 
1198.01 
1181.89 
1162.17 

^ 1145.90 

" 1133.77 
1127.68 

' 1111.93 



3- 13P(2) 1705.29 
2-llP(l) 1670.40 

1- IOR(O) 1655.40 

4- 12R(2) 1638.52 

2- 10R(0) 1621.05 

5- 12R(3) 1607.58 

3- 10R(0) 1588.98 
7-13P(l) 1579.18 



I- 8R(0 

0- 7P(l 
9 12P(] 

II- 14R 

1- 7R(0 
0-6R(3 
5-8P(l 

2- 6R(0 
0-5R(0 
10-9P( 
4-6R(2 
7-7R(l 
0-4P(l 

3- 5P(l 
9-7P(l 

0- 3R(3 

1- 3P(1 

1- 2P(1 

2- 2R(2 

3- 2R(l 
1-lPfl 



1544.94 
1518.90 
) 1503.78 
1) 1493.83 
1486.53 
1458.14 
1438.02 
1401.87 
1393.73 
) 1375.98 
1356.85 
1343.60 
1335.87 
1325.06 
1307.58 
1278.74 
1255.68 
1200.08 
1181.93 
1163.37 
1146.16 



15-5R(1) 1136.12 
2-lR(2) 1129.74 

3 1R(3) lllG.GG 



"Wavelength of tramsmissioii minimum at resolution R = 350. 
''AH transitions are Lyman band unless otherwise indicated. 
'^Werner band transition. 
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where we have taken Lgto = 2.5 x 10^^ ergs and (3 = —0.5. Therefore even if the GRB occurs in a 
low-density region, absorption by vibrationally-excited H2 may be observable if there is a molecular 
cloud on the line-of-sight within ~ 100 pc of the GRB. 

Early detection and accurate positions for GRBs will enable spectroscopic observations to be 
carried out while the fireball/afterglow is still bright. The High Energy Transient Explorer (HETE) 
is expected to provide rapid positions so that early ground-based observations may be possible for 
some GRBs. This may make allow ground-based telescopes to obtain R > 100 spectra of the 
afterglow for afterglows at z > 1.25 where the A = 1600 A region is redshifted into the U band. For 
rcdshift 2.2 < z < 4.3, the entire 1110 — 1705 A spectrum would be observable from the ground 
with silicon CCD detectors. 

The Swift Gamma Ray Burst Explorer mission, scheduled to be launched in 2003, is expected 
to return positions within seconds for ~1000 GRBs over its 3 year mission. The UltraViolet and 
Optical Telescope on Swift will image ~300 GRBs within 70 sec of the burst, with 6500-1700 A 
wavelength coverage. For the brightest bursts, observations through a grism can be made, yielding 
a 200-400 spectrum over the 6500-1700 A spectral range (Roming 2001). The UV sensitivity 
means that rest wavelengths 1705 - 1600 A (a range which includes 6 of the absorption features 
in Table 1) could be observed for GRBs at redshifts z > 0.06. For GRBs at 0.53 < z < 2.81, 
the UVOT spectral range will include the entire 1705-1110 A range where absorptions lines from 
vibrationally-excited H2 are conspicuous. 



The luminous optical transient associated with at least some GRBs will have drastic effects on 
surrounding molecular gas and dust, if the GRB is situated within a molecular cloud. To illustrate 
the effects, we have adopted a model light curve for the optical transient, and have considered 
power-law spectra Li, oc and Li, oc z/~^. We have calculated the resulting photoioinization 

and photodissociation of the surrounding gas, and thermal sublimation of the dust, for two values 
of the gas density, nn = 10^ and 10^ cm^^. Our principal results are as follows: 

1. Dust will be destroyed out to distances Rd ~ 10^^ cm in clouds of density nn < 10^cm~^, 
confirming earlier estimates by WDOO. 

2. The ionization and dissociation fronts are merged for spectra L^, oc and oc u^^. The 
destruction of the H2 is therefore due in part to photodissociation, and in part to photoion- 



8. Summary 



ization. 
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3. As proposed by Draine (2000), a substantial column density of vibrationally-excited H2 is 
created as long as the ionization-dissociation front has not broken out of the molecular cloud. 
This vibrationally-excited H2 will produce strong absorption lines in the 1705 — 1110 A region 
where normally Lyman-a would be the only strong absorption line. 

4. Spectra of GRB optical transients and afterglows should be compared to the predicted absorp- 
tion spectra shown in Figs. 9 - 11. Even if the signal-to-noise ratio does not permit detection 
of individual features, cross-correlation may make it possible to detect the ensemble of lines 
between 1650 and IIIOA. Detection of such absorption by vibrationally-excited H2 would be 
unequivocal evidence for association of the GRB with molecular gas. The UVOT grism on 
the Swift Gamma Ray Burst Explorer mission may be able to measure these features in the 
spectra of the brighter GRBs at redshift 2; > 0.1. 



We thank E. Roueff for making available H2 data, S. Hunsberger for providing the specifications 
for the Swift UVOT instrument, xxx for helpful discussions, and R. H. Lupton for availability of 
the SM software package. This work was supported in part by NSF grant AST-9988126. 
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